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1. CHARACTERISTIC DIFFERENTIAL EQUATIONS 


The differentiae, equations for the rotational ly symmetric potential 
flo-ws read 



If one starts from the velocity potential cp, and 

C 1 - - 2 + C 1 - - t = 0 (2) 

If one uses the stream function (see chapter I, equations (21) 
and (19) ) • Therein the velocity component In z direction is 


- — 

u = cp z = (1 - q. 2 ) k-1 I * 

r r 


(3a) 


the velocity component in r-direction 


v = cp r = - (1 - 






(3b) 


*'feotatIonsBymmetrische Potentialstromungen " (Eapitel IV ) , 
Technische Hoschschule Dresden, Archiv Nr. 44/4 November 18, 1940. 
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■while for the sonic velocity a in the selected non-dimensional 
representation 


a 2 = - q 2 ) (t) 

applies. q vas the amount of the velocity, thus '/u 2 + v 2 . 

In setting up the characteristic differential equations one starts 
heat from (l), since the derivation then is some -what facilitated- Later 
on, it is true, interest will he concentrated on the stream 
lines Tjr = const-, not on the potential lines cp = const j however, for 
potential flows the quantities cp or i|r enter only into the last 
characteristic differential equations whereas those treated first are 
completely free of cp or ijr and contain, aside from the independent 
variables, only u and v- u and v are more simply expressed by cp 
than by ijr . 

The characteristic condition (of. chapter II, equation (8), NACA 
TM 1242) is 


(l - + 2uvzr + - 2^z2 = 0 (5) 

If one now introduces the angle $ of the stream line and the symmetry 
axis (z-axis) of the flow, one obtains 


u = q cos -a v = q sin 9 (6) 

Due to the rotational symmetry, 9 may be assumed between 0 and k/2, 
if reverse flowB are excluded. 

Then the characteristic condition turns into 

(a2 - q 2 cos^iS)r^ + 2q^sin 9 cob aif + (a^ - q^sin^i3)z^ = 0 (j) 


f 

X 
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Two roots X ' and X ' ’ result from this relation for 4^, the slope 

dz 

of the characteristic hase curves with respect to the z-axis : 


q 2 sin 3 cos 3 + a\/q 2 - a 2 

2 2 2 
q cos 3 - a 


X 


q 2 sin 3 cos 3 


2 2 
q cos 3 



(8a) 
* (8b) 


As dif ferential equation for the first family of characteristic hase 
curves there results 


dv - X ' dz = 0 


(9a) 


for the second family 


dv - X ' ' dz = 0 (9*0 

(see chapter II, equations (28a), (29a), NACA TM 1242). 

For interpretation of these relations (9 a ) and (9b) - rather 
complicated due to the form, of X ’ and X 1 ’ - it is test to start- from 
the equivalent statement (7 ) • For the latter one may -write 


q 2 (cos 2 3 r^ - 2 sin 3 cos 3 zr + sin 2 3 z 2 ) = a 2 (r 2 + fc 2 ) 
or, with introduction of differentials, 

2 2 

(cos -0 dv - sin 3 dz) _ a_ 

. dv 2 + dz 2 q 2 


( 10 ) 


Then there is on the left side the square of the vector product of the 

unit vectors (cos -3, sin 3) and ( -=======, -======), which are 

Vvdv 2 + dz 2 Vdv 2 + dz 2 / 

tangent to the stream line or the characteristic base curve, respectively. 
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If one denotes the angle between stream line and characteristic base 
curve by a, (10) therefore states that 


sin 2 a = 4 (U) 

q 2 

This, however, is nothing else but the well-known Mach relation: a is 

the Mach angle between stream lines and Mach waves which are accordingly 
recognized as being identical with the characteristic base curves. 

Hence the interpretation of the equations (9a) and (9b) results 
directly. The positive z direction is assumed to be in the direction 
of the mean flow. Noting the signs of the roots in the expressions 
for X ' and X 1 ' one then finds that the first family of characteristic 

base curves (9a) foims with the stream lines the Mach angle a - arc sin 

to the left, looking in flow direct ion j the second family of character- 
istic base curves forms the same angle to the right. 



The first family of characteristic base curves in our specification may 
be designated as left-hand Mach waves, the second family as rirht-hand 
Mach waves . The second characteristic equation of the first family of 
characteristics in the form (chapter H, 28 b, MCA TM 1242) reads in our 
case : 

(l - ' du + (l - Z^dv + Idr = 0 (12) 

If one takes into consideration that In the expression (8a) for X ' the 

denominator is q 2 aoa 2 'd - a 2 = u 2 - a 2 , there results from (12) after 
multiplication by a^ and consideration of ( 6 ): 


-(q. 2 sin d cos d + a\Jq 2 - a 2 ^cos d dq. - sin dq. dd) ■ 

+ (a 2 - q^sin^d^sin d dq. + cos dq. dd^ + a^q. sin -d — = 0 


aj| a* 
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or, after the multiplication of the brackets has been carried, oat 
terms have been rearranged, 

|(a^ - q 2 )sin -3 - a\|q^ - cos aj'dq + ja^cos a + a|q^ - sin aj 


+ a^q sin 3 — = 0 


If one factors q from the first two terms, one obtains 
q(a 2 cos -fl + a\jq 2 - a 2 sin a^)(da - -S 


— - — dq ) + a 2 q sin a — 
aq ' 


= C 


or 


(cos a + ii 2 2 &2 sin a) (da - dq) + sin a = 0 


If one now presupposes that 


cos 


, ipr: » 2 

a + i 


S= 1 a sin a i 0 

a 


one may divide the equation (13) by thiB expression and obtain 

I 


da - fe-.. 8 - 2 dq + 
aq 


cot a + 


= 0 

'/q 2 - a 2 r 


Using the Mach angle, one may equate in (l*j) 


and the 
q da 


(13) 


( 1*0 


( 15 ) 
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Accordingly, one has for the second family of characteristics 


dd + \lz 2 ~ a2 

aq. 


dq + 


cot 



(17) 


with the presupposition 


cos -3 - ~ a - sin -3^0 (l8) 

What is the significance of the cases excluded "by (l^) and (l8), 
respectively? If, for instance, (18) does not apply, one obtains ■d = a. 
Since the stream line forms with the z-axis the angle ■a, with the 
characteristic base curve the angle a, the respective characteristic 
base curve then becomes parallel to the z-axis; thus one has on it at 
this point dr = 0, the last term in (17) would assume the indefinite 
form 00 x 0. 

One may avoid these difficulties completely by using instead of 
the differential dr in (15) and (17) the differential ds, the arc 
element of the characteristic base curves. Using dz, which suggests 
itself would, in contrast, only shift the difficulty elsewhere. 

If in the differential equation (15) perta inin g to the first family 

of characteristics the factor ~ is designated by E’. it is therefore 

r 

to be equated 


E 1 — = S’ M (19) 

r r 

with precisely the expression S’ to be determined. If one denotes the 
angle of the first family of characteristic base curves and the z-axis 
by P there is dr = sin p 1 ds. On the other hand, the slope of this 
characteristic base curve was given by X,' according to (8a), 

thus tan P * = X*. Therewith dr becomes dr = ds: thus, since 

there shall be B 1 dr = S 1 ds. 
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After a few transformations one obtains 


AH'* 


a sin -8 


sin d - \Jq^ - a 2 cos d 


)J\ + i _ q a sin d + \/q 2 - cos d 


q_ cos c d - a 


and finally 


S ' = — sin d = sin a sin -3 

a 


Correspondingly 


( 21 ) 


X"E" 


A s A fl -i. 


a sin d + ' ]j q5 - a^ cos 


d 



-1 


a sin -6 - \/q2 - a2 cos d 

2 2 a 2 

q. cos d - a 


The signs are selected so that, in agreement vith the geometrical 
derivation above, — ■ ' — = sin p * 1 (p ' ’ angle of the second family 

' ,2 + 1 

of characteristic base curves with the z-amis) . Therewith S' 1 finally 
becomes 


S’ ’ = - JL sin d = - sin a sin d ( 22 ) 

Thus one obtains without any distinctions of caseB as second 
characteristic equation for the first family: 

dd - cot a + sin d sin a — = 0 (23) 

<1 r 

and for the second family: 

dd + cot a - sin •A sin a ~ = 0 

q. r 


t 


( 24 ) 



8 


NACA TM 1244 


Only this simplification of the characteristic differential equations 
enabled us to perform the calculation of rotationally symmetric 
potential flows with tolerable calculation expenditure . To our knowledge 
no indication of the possibility of such a simplification is to be found 
in the suggestions of other authors (Fran k l, Ferrari) concerning bhe 
calculation of rotationally symmetric potential flows-*-* 


2. TREATMENT OF PRACTICAL EXAMPLES 


For the calculation of practical examples, two methods are at 
disposal for the approximate solution of the characteristic differential 
equation system, as was discussed in more detail in chapter U, 7 , 

NACA TM 1242 j we decided for the "field method" (chapter II, p. 16, NACA 
TM 1242, p. 16) . If one progresses from already known fields to new 
ones, the equations (cf. chapter II, formula (49), p* 23, TM 1242, p. 24 2 

A I X l'( p III " p l) + C Il(' 1 III “ 4 l) + D l( y III " ?i) = 0 ( 2 5a) 


■ A H X H "(^II " p Il) + C Il(^IH ~ *El) + " y u) = 0 ( 2 5 b ) 

are used, which result from the following characteristic differential 
equations (II, (28), ( 29 ), p. 14, TM 1242, p. 15 ) : 

AX - 1 dp + C dq + D dy = 0 (26a) 

AX"dp + Cdq + Ddy = 0 (26b) 

with X 1 signifying the slope of the first family of characteristic base 
curves, X 1 ’ signifying the slope of the second family of characteristic 
base curves ; the differentials are approximately replaced by finite 
differences (and differences of second and higher order then neglected) 


As noted afterwards, for once the arc element of the character- 
istic base curves has been introduced into the characteristic differ- 
ential equations from a purely mathematical side: (E* Eamke 's 

'iDifferentialgleichungen reeller Funktionen, " (Differential equations of 
real functions), Leipzig 1930, p. 389 (9))* We reserve for a later time 
looking into a possible connection of Karake 's expressions of the 
characteristic differential equations with ours. 

O 

mFrom now on, these references will be abbreviated; for this 
reference, for instance, we shall write: U (49), p. 23 , TM 1242, p. 24. 
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as has been described in detail before (II, p. 16 ff . and p. 22 ff . etc. 
'EM 12H2 p. l6 ff . and p. 23 ff .) I and II signify the numbers of the 
fields with already known approximate values of p and q whereas III 
denotes the adjoining field for which these approximate values are to 
be determined. 

The relations (25), derived from the quite general form of a 
quasilinear differential equation of the second order (II, l(l), p. 2, 

TM 1242, p. 1) are now to be given for a special use in adjustment to 
the flow type of importance . 

For the rotationally symmetric potential flow there correspond to 
the two equation (26) according to (23) and (24) the pair of equations 
(second characteristic equation for the first and for the Becond family 
of base curves) 


d-fl - cot a + sin 9 sin a — = 0 
1 r 


(27a) 


d-3 + cot a - sin a sin a — = 0 
9. r 

in addition, the Mach relation (first characteristic equation) 


(27b) 


sin a = — 

a 


( 28 ) 


is used. Therein a was the local sonic velocity, q the amount of 
velocity, a the angle of the stream line and the symmetry axis (z-axis) 
a. the Mach angle between stream lines and characteristic base curves 3 , 
r the distance of a point from the z -axis . 


Within the scope of our approximation principle we write for (27a) 
and (27b) in analogy to (25a) and (25b) - putting for abbreviation 


cot <x v \ q^ - a^ 




V^V 


= Q V 


(29) 


3 Cf . remark on p. 4 
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and 


sin sin Oy = Sy (v = I, II) 


(30) 


" *i " q i(^iii ~ ^ + s i 


s iii " s i 


= o 


(31a) 


(progressing along a left-hand family vith the angle $ + a toward 
the z axis) and 


%n “ %i + %i(ta " <iii) " " 0 ( 31b ) 

respectively (progressing along a right-hand family with the 
angle - a toward the z - axis); Sjjj - s v simply signifies the 

distance of the field centers IH and V from each other. 


Following the formulas required are given; they are to he used 
in connection with the schemes given hy figure 1.. 


From the two linear equations (31) q.jjj (and therewith, according 
to (28), also Ojjj) and dj-j-j for the unknown field III are easily 
found. The el im ination of t3 tt t from (31a) and (31h) first yields 


Qill + Qxiln + s i 


*111 


B III - B I 


Qj + Qji 


+. s. 


II 


B i n - sn 
r ll 


+ S U - *1 


(32) 


hy insertion of Ijjj into (31a) or (31b) then results 

fl m = *1 + Q l(^-III " %l) " ®i — - 81 (33a) 

or 

*111 - *n - Q Il(^III - ^II) + s n — ^33b) 

Of the two formulas (33) one is, of course, superfluous; however. It may 
he useful for control of the calculation. 
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The relations ( 32 ) and (33) are applied in the standard case 
(cf • fig. la), where one deals with Inner fields only; if, on the other 
hand, one reaches, in progressing, the fixed "boundary (case of fig. lb 
and lc, respectively), -where is prescribed by the boundary condi- 
tion ^), one has, according to ( 31 ), the simpler formulas 


a ni 


s ni - St 

Qjqj + s i rx 1 + ~ *i 

0,1 


(case lb) (3^a) 


and 


iin ■ 


Q n q -n + s n 


s iii - s n 

r n 

S:i 


+ ■a 


11 


hi 


(case lc) (3^b) 


respectively . 

Conditions at the boundary of a free jet are quite analogous. 

Only is now not prescribed but, on the contrary, to be determined 
■whereas pressure and therewith the amount of velocity q are known. 

The solution of (31a) and (31b) with respect to 3 yields^, for given q 


*111 = *1 + Q i(lni 

and 

*111 = *11 ' Q n( q -m 

respectively 


“ *l) - ^ 


S HI " S I 


(case lb) (35a) 


- q.jj^ + S^-j. -SI (case lc) (35b) 


The formulas so far refer - as will generally be the case - to 
progressing in flow direction. Now and then, however, (cf . our procedure 
in the second example, Bection It) it is also required to follow the flow 
rearward and one then has the conditions sketched in figure 2. 


\Ef necessary, one assumes the boundary approximated by series of 
straight lines 3 = const. 

A peculiarity here experienced in the determination of the field 
center will be treated in detail in the discussion of our second 
example . 
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The corresponding formulas naturally follow directly from the 
former ones, one only has to replace the progressing from field I to 
field III "by that from IH ' to II ' and, in analogy, the progressing 
from H to HI by that from III ' to I ' . Solution of the new 
equations with respect to Ijjj ■ and ^xil 1 results^ in 


q -Hl' = 


Qjilj.+Qjjrqjjt -Sj, 


3 III' 


r r 

Q, ^ i + Q' 


-S H 


3 II 


i - S' 


HI ’ 


L II 


' a n ' + 


II 1 


and 


(36) 


J m' “ fl i* + Q i ' " ^ii') " s i' 


B I* “ a IH 


r T » 


(37a) 


a HI' “ *ir' " ^II'^H' " %Il) + 


B-__- , - a 


rr 1 in ' 


II ' rjj , 


(37b) 


respectively) 


(case 2a) ) 

S T 1 “ S TTT 1 „ . 

Qi.q]-, - Sj, — — *ni' + * x , 

n , s ■ 1 ^ - — ■ ■ — ■ - ? d , = d 

HII ' n ) HI ' boundary 

i 

( 38 a) 

(case 2b )j 
^I ' " S HI 1 

Q II '^II* " S H' rn> fl II' + *111 ' 

^III' = Qjj, 31111 = ^boundary 

(38b) 

(case 2c) . 


For the free jet formulas for the progressing upstream were not 
required and therefore not derived. 
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Furthermore, knowledge of the following two additional remarks is 
important for the application of the formulas above • 

1. For the rotationally symmetric flow it is sufficient to construct 
for a longitudinal section the flow pattern above the z-axis -9=0 
(■which may be Interpreted as boundary), since the, course below this 

axis is simply obtained by reflection on it. However, if one wants to 
calculate this course separately - for instance for purposes of control -, 
one must take into consideration that one is now dealing with a left- 
hand Cartesian z-r system. This is taken into account by providing in 
the formulas the quantities r v - in themselves, according to definition, 
positive - with the negative sign. 

2 . If one uses in the formulas a field center V lying on the 
symmetry axis -9=0, the're becomes for it Sy = 0 and simultane- 

Br, - S „ 

ously r v = 0 ; the expression Sy-^ then assumes an indefinite form. 

r V 

In this case the significant expression is simply to be omitted, since 

there results lim — = 0, for the following reasons: Along the stream 

r — * 0 r 

line d =0 there is v = 0 and, for reasons of symmetry, obviously 

also — = 0 . Thus one has, since v = q sin -9: 
dr 

0 = £9. sin 9 + q cos 9 — and consequently — = 0 for r = 0 . Now 
dr dr dr 

one may "write for small values A9 and Ar : — = P" x B AP ^ ft sin a — . 

r r Ar 

For transition to the limit Ar — *0 there results, because of = 0 , 

dr 

also ^IItHq — = 0 . 

Moreover, the following has to be said in general for practical 
application of our method. It will be useful to represent, before 
starting the calculation proper, the values of Q, and a as functions 

of q by means of the relations (4) a 2 = K r ... jr (l - q 2 ^ and ( 29 ) for " 

the interval a 2 = . K 3- <_q ill: of course, this has to be done for 

Y k + 1 

every a value separately. For this purpose a table was set up which 
indicates the Q, values and a values for q steps of 0 . 02 ; for the 
intermediate values linear interpolation which is easily performed is 
sufficient. It was shown for the used K-values (k = 1.405; k = 1.18) 
that for the chiefly important q domain (if q does not approach the 
unit too closely) Q may be practically regarded as linear function 
of q. The two tables for k = 1.405 and k = 1.18 on which the 
calculation of our examples is based are added at the end (pp. 26 and 27 ). 
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The quantity S, on the other hand., is determined from case to 
case. Its calculation may "be performed very simply "by slide rule 

according to (30), all the more so since the factor sin a = — ■which. 

a 

appears as intermediate quantity in the calculation required for the 
setting up of the table may be taken directly from that table • 

For the rest, it will be noted in using our formulas that certain 
combinations of quantities occur repeatedly, which means a saving in 
calculation expenditure. 

The calculation and construction of the field HI - which for us 
is to be representative for all fields nevly to be determined (cf . 
fig. 1) - is performed as follows: first, q. In and are deter- 

mined according to the respective formulas, with the quantities denoting 
length Sju - s v and r v to be taken from the drawing (its scale is 

of no importance since only the quotient of those lengths enters into 
the calculation) . The table yields the pertaining values of ftjjj 

and O q~pr , and with the aid of one can find S T TT . (it haB to 

be noted that Q, and S in themselves are auxiliary calculation 
quantities which are required only for the continuation of the proce- 
dure . ) There follows a geometrical process, the closing of field III 
which requires knowledge of the Mach angle ajjj and of the stream 
line angle a -nrr • At the comer point 1 (cf. fig. la) the 
angle + Oqjjj, at the corner point 2 the angle a jjj - ccj-q- is 

subtended by the horizontal ( | | z-axis) and the legs of both angles are 
made to intersect at the desired comer point 3 • Since thereby, 
naturally, a certain moment of inaccuracy (aside from the errors caused 
by the approximation principle) enters, the position of the corner 
point 3 "was occasionally checked analytically by means of the equations 
(II (42), p. 18 , TM 1242, p. 20) . 

r 3 “ r l = ( z 3 “ z i) x tan f a TTT + ccq-jj^ 

, X v (39) 

r 3 " r 2 - ( z 3 " z 2 ) x tan (%n - “nij 

or, solved with respect to the coordinates z~, r„: 

a * S 


Z 1 x + °m) - zgtan^jjj - ctj-jj) - r 1 + r 2 

tan (a in + ccq-jj 'j - tank’d jjj - 


( 39 ') 
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After this general representation of our combined calculative- 
graphical method we now turn to the description of the special examples 
that have been carried out. 


3 • FIRST EXAMPIE i Dili' USER 


This example "which we selected ourselves and which represents the 
first attempt of application of our method, deals with the flow of 
air (k = 1.405) through a conical nozzle of the half opening angle 6° 
according to figure 3 , which shows a longitudinal section. As initial 
condition we assumed q. in this longitudinal section on a circular arc 
(the apex of the cone is the center of the circle), to be as constant 
throughout, equal 0.59> (the pertaining Mach angle is 38° 2'); "3 results 
from the fact that the stream lines must intersect this circular arc 
radially. The distribution of the fields is here selected still 
relatively roughly (-a step of about 4-° in r direction) : field comers 

on the circular arc at 6=0 and 6 = ±4°, thus field centers 
at 6 = ±2° and 6 = ±6° (boundary), as can be seen from figure 3* 

Also, the method was extended* only to the construction of a few fields 
since no peculiarities appear. 

The construction in general has already been described in detail; 
supplementary information is needed only about the determination of the 
field centers. For inner fields one finds the new field center III 
(cf . fig. la) as center of the diagonal 12 of the comers 1 and 2 of the 
known fields I and H. At the boundary the procedure is as follows 
(cf. fig. lb and, analogously, fig. lc) : the parallel to the character- 

istic base curve 12 of the field I is made to intersect with the 
boundary line, and the intersection point is defined as field "center" 
of the half -field HI. 

The ratio of the prevailing pressure p to the tank pressure p Q 
(that is, the pressure for which the flow velocity equals zero) has been 
plotted into the individual fields . This ratio results, according to I, 
2, from the respective q. according to the equation 

K 

-B- = (1 - q 2 )"' 1 (40) 

Po 

Furthermore, a number of stream lines are drawn; they were selected 
according to the 'principle that in every stream tube the, same mass 
transport takes place per unit time • These stream lines are attained by 
progressing by equal amounts As]r, starting from 6 = 0. Their position 
is fixed in the following ma nn er. 
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The point which -was taken in figure 3 as starting point of the 
stream line 3 = 0 is assumed to he the origin of our z-r coordinate 
system. One now searches along z = 0 for the r or din ate of the 
points through which the stream lines must pass . For constant z one 
has, if the first equation (l6) of chapter I, p. 9 is written as 
difference equation, 

1 

. Aijr ft r q cos £(l - q?) K ^ Ar 

»P Ar (4l) 

According to the Simpson rule? 


P r=r boundary „ 

^boundary " ^o - ^ ^ ^ 


(42) 


was numerically determined (which requires setting up of a table for a 
number of r values and of the pertaining values of P) and 


A'if = 




boundary 

n 



(43) 


was then found, with n signifying the number of stream lines to be 
plotted (in the present case n = 6) . Insertion of Ai{r in (4l) then 
yields the n different Ac by which one must progress along z = 0, 
starting from r = 0. Xt is practical to proceed from the boundary: 
one introduces P(r = ^boundary) in- 4° (4l) and thus obtains Apr. Then 
one determines by interpolation from the table the P value pertaining 
to r bo undar y ' A l r 811(1 inserts it again in (4l), thus obtaining Agr. 

The same procedure is performed for P = P (boundary “ Apr “ Agr) an* 1 

n 

yields A^r, etc. Naturally one must not expect !EZ A^r to become now 

V=1 

exactly equal to r t omi aar y, sixLCe fl nst, always a too large P was 

used and, second, P waB regarded as constant for the entire (finite) 
interval Ar. Within the scope of graphical accuracy, however, it 
proved sufficient to plot now r against A^, to alter the -curve drawn 
through the ri (i = 1 , 2 , . . ., n) so that it passes exactly 
through r = 0 for n Aty = ^boundary " ^o (insignificant correction’), 
and to use the new coordination from the r-AJt diagram. Besides, one 
may also, as a proof, proceed inversely starting from r = 0; one then 
obtains too large values for Ar, whereas those found according to the 
first method were too small. 

7of course, another suitable quadrature formula could have been 
selected just as well. 
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The stream lines start at the angle -3 given by .the Initial 
condition (the slight difference between the circular arc and the 
rectilinear section z = 0 could he neglected) ; -when they reach the 
field boundary, they are continued at the new inclination coordinated 
to the adjoining field and thus result in the broken lines shorn in 
figure 3 . 


4. SECOND EXAMPLE: NOZZLE 


Here the flow of a gas with k = 1.18 through a nozze of pre- 
scribed longitudinal dimensions, the shape of ■which in longitudinal 
section is given in figure 4, was to be determined. The nozzle is a 
truncated circular cone with the wall angle 6°10 ' . The narrowest cross 
section has a diameter of 118 millimeters, the exit cross section, one 
of 178 millimeters. At the narrowest cross section the obtuse cone 
starts with a slight - not further defined - curvature . The location 
of transition is marked in our drawing. The tank pressure is 14 atmos- 
pheres, the pressure in the exit cross Bection 1 atmosphere. The 
determination of the curve results necessarily in our procedure which 
conforms to a report of Th. Meyer® and will be described in detail 
below. 

In the surroundings of the narrowest nozzle cross section where 
the flow velocity is near Critical velocity, our method is not yet 
usable in practice (and this applies to the supersonic region, also), 
because the Mach angle a still differs too little from 9°° and. the 
fields thus become much too narrow to allow even a tolerably rapid 
progress in the construction of the flow pattern. For this reason we 
apply our method only at a sufficient distance from the narrowest cross 
Bection (naturally making our arrangements so that we need operate only 
with the rectilinear boundary) j it is true that initial conditions 
for q. and S for the start then must be obtained by another method. 

This latter, accord in g to Th. Meyer, consists in assuming at the 
point of the z axis where sonic velocity is passed as origin of the 
coordinate system a series development with respect to the two 
coordinates z and r for the velocity potential 9 . It is pre- 
supposed that in the proximity of the origin the first series terms 
already give a good approximation; proof of this will be given later 
(cf. p. 23). The coefficients of this development result by insertion 
of the series into the differential equation for 9 (I ( 21 ), p. 12 ) 


^Meyer, Th.: Tiber zweidlmensionale Bewsgungsvorgange in elnen Gas, 

das mit ijberschaUgeschwindigkeit stromt. Diss. Gottingen 1907, 
published - in Mitteilungen uber Forsch.- Arb. ing.-Wes . Heft 62, 
Berlin 1908. 
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and eut sequent comparison of coefficients. However, certain coeffi- 
cients remain arbitrary as long as no further conditions are added. The 
first condition is obviously that the flow should be symmetrical to 
the z axis, that is, v must not contain other than odd powers of r 
or cp other than even powers of r. As second condition Meyer intro- 
duces the assumption (well confirmed by tests) of a linear velocity 
increase along the z axis (with the critical sonic velocity at the 
origin) 7 

A priori, another method of determining the coefficients would be 
more readily apparent: .it is the method used by G. I. Taylor. Taylor 9 
and with him S. G. Hooker 1 -*- 5 who transferred Taylor's method to the 
three -dimensional case of a nozzle of circular cross section, do not 
have Meyer 's second condition but prescribed for their calculations a 
certain nozzle with circular-arc boundary. The boundary conditions 
then yielded the lacking qualifying equations for the coefficients. 
However, it becomes evident that in this case the determination of the 
coefficients becomes more complicated than according to Meyer's method, 
quite aside from the fact that the latter method is better adapted to 
our statement of the problem where the form of the boundary is not 
accurately fixed. Moreover, the results of T. E. Stanton's-^- measure- 
ments on the nozzle selected by Hooker confirms the usefulness of 
Meyer's assumption (see above), that is, one obtains in this manner 
reasonable nozzle shapes. 

If one denotes all quantities made non-dimensional by a dash above 
them and refers the velocity to the critical sonic velocity a*, Meyer 's 
second condition reads: 


u = 1 + z 


(44) 


%aylor, G. I.: The Elov of Air at High Speeds past Curved 

Surfaces. Aeronautical Research Committee Reports and Memoranda 
Wo, 1381, London 1930. 

^Hooker, S. C.: The Flow of a Compressible Liquid in the Neigh- 

bourhood of the Throat of a Constriction in a Circular "Wind Channel. 
Proc. roy. Soc . London A, Bd. 135 (1932), pp. 498-511* 

Cf. also Gortler, H.: Zum Ubergang von Uhtershall=r zu TJber- 

schallgeschwindigkeiten In Eiisen. Zeitsch. f. angew. Math. Mech. 

Bd. 19, Nr. 6, Lez. 1939 • 

^Stanton, T. E.: Velocity in a Wind Channel Throat. Aero- 

nautical Research Reports and Memoranda No. 13 88, London 1931* 
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The stream lines satisfy the relation 

M = || = I (J+5) 

dz 1 cpz u 

with r expressed as power series in z, the coefficients of which are 
to he dete rmin ed according to (45) « 

If one proceeds up to terms of the fourth order, the expression 
for the velocity potential becomes - with the first condition taken 
into consideration 

cp = ajz + a 2 Z 2 + cgr 2 + a3z3 + b 3 zr 2 + + b 4 z 2 r 2 + c]^ + . . .; 

(46) 

because of the second condition one has a-j_ = 1, ag = i, a^ = a^ = 0, 
consequently (46) may be replaced by the simpler form 

cp = z + ^z 2 + CgT 2 + b^zr 2 + b^z 2 ? 2 + c^r4 + . . . (47) 

with the four unknown coefficients Cg, b^, b^, c 4*. Hence results 

t 

u = || = l+z + b^r 2 + 2b^zr 2 + . . . (48a) 


and 


v = = 2 c 2 X + 2b^zr + 2b^z 2 r + 40^3 + 

Insertion into the differential equation I (21) 


(a 2 - u 2 ) cp Eg - 2 uv$j; E + (a 2 - v^cp^ + £ cpf = 0 


(48b) 



with 
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leads by comparison of tbe coefficients to the four equations: 


+ 

H, 

K. 

11 

0 

(49-1) 

+ l) X 2 a* 2 = 0 

(49-2) 

(k + l)a*^ = 0 

(49-3) 


(k - + 1604 - l6a*^c^ - 4 c 2 2 J- — |g X 1602 ^ - (k + 1) 

X ^2a*^b^ + Sa*^b^ + 8 02 ^ 3 ) = 0 
They yield for the unknown coefficients the values: 

C2 = 0 


3 2(l - a" 2 ) 


b 4 = 1. + J a . 


*2 


4(l - a* 2 ) 2 


°4 = 


(49-4) 


(50.1) 

(50.2) 

(50-3) 

(50-4) 


read 


I6(l - a* 2 ) 2 

Let the equation of the stream line up to terms of fourth order 

F = a + Pz + 7 z 2 + 8 z 3 + ez^ (51) 

If one enters with (51) into (45), there follows according to (48) : 

o o 2cof + 2b -Zr + 2b\,z^r + 4ci,r3 

P + 27Z ->■ + 4 ez J = ^ - -jr" ■£— 

1 + z + b^i^ + 2b] + zr^ 


( 52 ) 
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If one inserts, moreover, on the right side of (52) the value 
of r from (51) , the comparison of the coefficients yields the 

folio-wing four relations^ 2 : 



P(1 + 

b 3 a?) = 

4c4a^ 

(53.1) 

27! 

^1 + b^a?) = 2b + 

(lSc k - 

2b^)a 2 3 - 0 - 2b^c c(3 2 

(53 -2) 

38^1 + b 3 a 2 ) 

= 2bj0 + 2bl^a - 2 7 

- * 3 f> 3 

- 6b 3 a£7 -(Vb^ - 22c^)x(aJ3 

2 + a 2 7) 





(53-3) 

4e(l + h 3 a 2 ) 

= 2b 3 7 + 2b - 3S 

- 8b 3 a£8 - 4b 3 ^p 2 7 + ay 2 ) 



- (l2b^ - 240^(0^7 + V + IT./ (53-4) 


The quantity a which indicates the r value (that is, half the 
nozzle width) at the point (z = 0) where at the nozzle center the 
critical sonic velocity a* prevails, remains, as an integration con- 
stant, at first arbitrary, a will then serve for the determination 
of the scale which is used to mate z and r non-dimensional. 

We may now turn to the filing of the curvilinear boundary section 
near the narrowest cross section of our nozzle (cf . p. 15) which had so 
far been left undetermined . This curve section must be stream line and 
must pass with continuous tangent into the rectilinear section with the 
prescribed inclination of 6°i0 ' . Furthermore it must be observed that 
agreement with the prescribed dimensions is maintained and also that 
the z -interval Is not exceeded (only within this z -interval the 
admissibility of our series expression may be regarded as satisfied) . 

By trial an a value was determined which Is coordinated to a stream 
line suitable for our purpose. 

The performance of the numerical calculation yields for tt = 1.18 
according to (50): 

b 3 = 0.545 hi,. = 0.3705 = 0.0743 


12 

At this opportunity, a group of misprints in Meyer's report 
should be pointed out. On p. 62 one must read in his corresponding 
equations (always first term on the left side): in (2), 27 instead 

of 7, in (3), 3S instead of 5; in (4), 4e instead of in (5), 
5£ instead of in (6), 6t) instead of tj. 
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For a we selected the value 0*5 (it should he noted that this 
value must not he compared with Meyer 's a = 0 .5, Meyer having made his 
calculations for two-dimensional air flows with k = 1.405, hut that 
its corresponding value in Meyer's report - with respect to the course 
of the stream line - is prohahly about a = 0.2). Then (53) yields: 

P = 0.03265 7 = 0.2258 5 = -0.02677 = 0.0651 

The equation of the boundary section ( = stream line) reads 
therefore, according to (51) : 

r = 0.5 + O.03265Z + 0.2258z 2 - 0.02677z 3 + 0.0651^ (51*) 

From the condition 

= tan 6°10 1 = 0.1081 = 0.03265 + 0.45l6z - 0.08031z 2 + 0.2604z"3 
az 

results for the z at the Junction (in fig. 4 marked on the boundary) 

z q — 0 . 1695 

To this z 0 pertains, according to (51*), the value 

r Q = 0.5119 

ct = 0*5 corresponds to the radius at the entrance of our 
nozzle r Q = 59 millimeters j thus ve must multiply the above values of z 

and r by = 118, in order to find from Meyer's quantities the true 
0-5 

lengths in mini meters for our example. 

If one inserts the numerical values of the coefficients in (47) 
and (48), one obtains (a* = O.287) 

<P = — = z + |^ 2 + 0.545zr 2 + 0.3705z 2 r 2 + 0.0743r^ (47*) 

u = —gg- = 1 + z + 0 •545r 2 + 0 •74l0zr 2 (48a*) 

v = = 1.090zr + 0.74l0z 2 r + 0.297f3 


(48b*) 
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With the aid of the equations (48*) one could finally, according 
to the relations ‘ 


q. = \Ju2 + v2 $ = arc tan — 

’ u 

calculate along the rectilinear section z = z Q the initial values 
of q. and -o "which "were requirements for our method. 

The construction of the flow pattern in the interior of the nozzle 
to its exit proceeds in analogy to the construction of the first example 
described before in detail; only the subdivision of the fields -was now 
selected finer than in that example . Finally it should be mentioned 
that, starting from z = z 0 , we constructed the flow pattern also for 
some length rearward and found the q. values determined according to 
our method in good agreement with the corresponding values resulting 
from the Beries development . 

Figure 5 shows the gas flow after leaving the nozzle for the case 
of an external pressure of 1 atmosphere. Figure 5 is the direct 
continuation of figure 4 and to be considered joined to it at the right. 
Since it could be seen that the subdivision of the fields in the nozzle 
interior was no longer fine enough for the construction of the free jet, 
it was made from new initial conditions for q. and -a • These values 
were obtained by interpolation from the course of q. and -a (discon- 
tinuous due to" our approximation method) along the straight 
line z = constant at the end of the nozzle. 

The behavior of the Jet depends on whether it enters, after 
leaving the nozzle, a region of higher or lower pressure. In our case 
the pressure at the boundary of the nozzle exit was found to be larger 
than the external pressure; it is true that the pressure difference 
proved to be only slight. Figure 5 indicates that P a / p Q = 0.085; 

since p Q = 14 atmospheres, there follows p a = 1.2 atmospheres. Thus 

there results as difference: 1.2 atmospheres - 1 atmosphere =0.2 atmos 

pheres. The boundary point forms the starting point of a rarefaction 
wedge with the opening angle ct^ - <x a (with signifying the Mach 

angle pertaining to the internal pressure, the Mach angle pertaining 

to the external pressure) which extends into the interior of the Jet. 
Here this angle was so snial 1, (an ample half degree) that we could, 
within the scope of graphical accuracy, regard the rarefaction wedge as 
practically infinitely thin and coinciding with the characteristic base 
curve starting at that "singular" point. Still, one can see clearly in 
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figure 5 that the stream lines along this curve experience *a break 
exceeding the customary amount caused by the approximation"^. 

Our construction deviates from the foregoing due to the fact that 
we deal no longer with a fixed boundary but with a free boundary to be 
determined, with the q^ pertaining according to (it-O) to the external 
pressure prescribed as boundary condition. First, the continuation of 
the straight boundary line (d = 6°10 ') at the singular boundary point 
has to be determined, which is done in principle by formula (35a) * 
However, both field centers I and HI now coincide, precisely at 
the boundary point; consequently ds becomes zero, and the third term 
of the sum at right is eliminated. Since we may, as was mentioned 
before, assume Q in the interval in question q.i . • . q. a as linear 
function of q, -fl a may be found from <8* = 6°10 1 without step-by-step 
difference calculation according to (35aJ dir ectly by an integration 
process : 


= fl i + 



Q dq. 


(52) 


with Q, = aq + p(a, 0 constants which can easily be determined from 
the established table (cf. p. 26 )). 

If one wants to proceed further along the boundary, a new principle 
must be introduced for the geometrical determination of the field center 
on the boundary, since the one indicated for the fixed boundary now 
fails to work. Let the field I be known (of. fig. 6) . Point III is 
then fixed by elongating the line that so far formed the boundary (the 
first time the direction is given by -s a calculated according to (52)) 
beyond the comer point I and having it intersect with the parallel to 
the section 23 * One now lets the new direction t9ju found by calcu- 
lation from (35s) start at III. The boundary stream line of the free 
Jet therefore undergoes the ( discontinuous) variation in its direction 
in the field center, in contrast to the inner stream lines where this 
variation takes place along the field boundaries . 

We had actually intended to pursue the course of the Jet until we 
had exhibited the periodicity in the expansion-contraction process. 
However, having progressed (on the outside) by about the length of the 
nozzle, we were forced to break off our procedure by an incident the 

— -jo ' ■ 

This deflection in itself takes place (in continuous transition) 
in the interior of the wedge and must be constructed (subdivision 
of - a a into a number of Ax-steps) . 
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interpretation of -which shall he postponed for the moment* It -was found 
that t-wo adjacent characteristics intersected (the point of intersection 
is outside of the range of fig* 5j however, the two curves in question 
are easily recognized) * The problem arises whether we have to deal here 
only with an apparent singularity - caused by the inaccuracy of the • 
approximation - or whether actually a compression shock occurs, although 
only a weak one* The latter is probable since we constructed the sec- 
tion of the free jet given by figure 5 twice and found in both cases an 
intersection of two adjacent characteristics moreover, on the basiB of 
our tests made in this direction one cannot reject the supposition that 
the same phenomenon will repeat itself at other points later* 

As a conclusion it should be mentioned with regard to the second 
example that the ratio p/p G again, as in the first example, was 
written into the individual fields (cf . p. 15) and that the 24 stream 
lines have been selected according to the same principle. 


Translation by 
Mary L* Mahler 
National Advisory Committee 
for Aeronautics 
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APEENDH 


TABLE I 

SONIC TELOCITY a, % AUXILIARY QUANTITY OF CALCULATION Q 
Q 

AND MACH ANGLE o AS FUNCTIONS OF TEE VELOCITY AMOUNT q. 


FOR k = 1.405 


4 

a 

a 

4 

0 A 

11 - (if 


a 

0.411 

0.4328 

1.000 

0.000 

90 0 

.43 

.4062 

.946 


.806 

71°5> 

.45 

.4019 

.893 

1.320 

63°15 ’ 

.47 

•3972 

.845 

1.344 

57°37 ' 

•49 

•3920 

.800 

1.529 

53°5 ' 

•51 

.3870 

•758 

1.682 

4 9 ° 17 i 

•53 

.3813 

.720 

1.816 

. 46°2 * 

•55 

• 3757 

.683 

1.946 

4 3 °5« 

•57 ' 

•3693 

.648 

2.061 

4o°24 ' 

•59 

.3630 

.616 

2.172 

38°2* 

.61 

•3562 

.584 

2.275 

35044 . 

•63 

.3498 

•556 

2.379 

33047 , 

•65 

•3415 

•525 

2.490 

31°40 1 

•67 

• 3338 

.4980 

2.599 

29°52 ' 

•69 

.3256 

.4720 

2.710 

28°30 ' 

•71 

.3170 

.4466 

2.824 

26°32 ' 

•73 

.3076 

.4213 

2.946 

24°56 ' 

•75 

•2977 

.3969 

3-080 

23°22 1 

•77 

.2870 

• 3730 

3.231 

21°53 ' 

•79 

.2760 

.3495 

3.395 

20°28 ' 

.81 

•2639 

.3257 

3-586 

19° 

•83 

.2510 

•3023 

3-799 

17°36 ' 

•85 

.2370 

.2789 

4.051 

l6°12 ’ 

•87 

.2220 

•2552 

4.352 

14°48 • 

•89 

.2051 

.2306 

4.740 

13°20 ' 

•91 

.1866 

.2050 

5.240 

n°5o l 

•93 

• 1 655 

.1780 

5.945 

10°15 ' 

•95 

.l4o4 

.1478 

7-04 

8°30 1 

•97 

.3096 

.1130 

9.08 

6°30 ’ 

•99 

.0638 

.0644 

15.67 

3°42 ' 

1.00 

.0000 

.0000 


00 

0° 
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TABLE H 

SONIC TELOCITY a, p AUXILIARY QUANTITY OF CAICULATION Q 
AMD MACH ANGLE a AS FUNCTIONS OF TEE TELOCITY AMOUNT q. 


FOR K = l.lfl 


a 

a 

a 

q. 


a = arc sin 

it - 

r> _ V a 

a 

0.287 

0.2874 

1.0000 

0.000 

90° 

.30 

.2862 

.954 

1.048 

72030 ' 

•32 

.2843 

.889 

1.610 

62°40 ' 

•34 

.2821 

.830 

1-977 

56°5* 

•36 

.2796 

•777 

2.251 


.38 

.2775 

.730 

2.462 

46°55 1 

.40 

,2750 

.688 

2.640 

43°25' 

.42 

.2724 

.649 

2-793 

40°25' 

.44 

.2693 

.612 

2.936 

3704a 1 

.46 

.2663 

•5795 

3.058 

35028’ 

.48 

•2633 

.5490 

3.172 

33°17 1 

•50 

.2598 

•5195 

3.288 

31°20 ' 

•52 

. .2563 

.4932 

3.397 

29033 . 

•54 

•2530 

.4686 

3.491 

27°56 ' 

•56 

.2406 

.4440 

3.608 

26°21' 

•58 

.2445 

.4218 

3-710 

24057 » 

.60 

.2400 

.4000 

3.819 

23034* 

.62 

.2354 

.3800 

3.931 

22°20 ' 

•64 

.2304 

.3602 

4.050 

2l°12 ’ 

.66 

.2254 

•3417 

4.168 

19°58' 

.68 

.2200 

.3236 

4.300 

18°53 ’ 

.70 

.2142 

.3061 

4.447 

17°50’ 

.72 

.2081 

.2891 

4.601 

16°49 ' 

•74 

.2017 

.2722 

4.770 

15°48 ’ 

.76 

•1949 

.2563 

4.961 

14°52 ' 

.78 

•3877 

.2406 

5-17 

13055, 

.80 

.1800 

.2250 

5.41 

13° 

.82 

•1717 

.2093 

■ 5.70 

I2°5 1 

.84 

.1627 

•1938 

6.02 

U°10' 

•86 j 

• 1531 

.1781 

6.42 

10°15 1 

.88 j 

• 1425 

.1620 

6.92 

9°20 ' 

.90 ! 

.1308 

.1454 

7.56 

8°22 ' 

•92 ! 

.1176 

.1278 

8.43 

7°20 ' 

-94 

.1023 

.1088 

9.72 

6°15 ’ 

•96 

.0841 

.0876 

11.84 

4°50 ’ 

.98 

i 

•0597 

.0609 

16.74 

3°26’ 

1.00 

.0000 

.0000 

ca 

0° 


til cH 



TABLE HI 


EEESSORE p HEffEBEED TO r nTBi EAHK 
PRESSURE p 0 AS A IOTCTION CCF 
THE VELOCITY AMOUNT q. 


4 

P/Po 

k = . - 1.405 

p/p 0 

K » 1.18 

0.30 


0.5388 

.32 


.4924 

•34 


.4467 

.36 


.4024 

•38 


•3597 

.40 


•3190 

.42 

0.5100 

.2800 

. 44 - 

.4741 

.2440 

.46 

.4382 

.2102 

.48 

.4030 

.1796 

.50 

.3686 

.1518 

•52 

•3350 

.1265 

.54 

.3023 

.1018 

•56 

.2710 

.0849 

•58 

.2410 

.0682 

.60 

.2127 

.0536 

.62 

.1859 

• 04 l 6 

.64 

.1608 

.0316 

.66 

•1373 

.0235 

.68 

.1161 

.0171 

•70 

.0966 

.0121 

.72 

•0791 

.0083 

.74 

.0638 

.0055 

.76 

.0502 

.0035 

.78 

.0386 

.0021 

.80 

.0289 

.0012 

.82 

.0207 

.0007 

.84 

.0144 

.0003 

•86 

.0094 

.0002 

.88 

.0057 

.0001 

.90 

.0032 

.0000 

.92 

.0015 

.0000 

•94 

.000 6 

.0000 

.96 

.0002 

.0000 

.98 

.0000 

.0000 

1.00 

.0000 

.0000 







